Effects of the Stark shift on the evolution of the field entropy 
and entanglement in the two-photon Jaynes- Cummings model 


Mao Fa Fang 

Department of Physics, Hunan Normal University, 
Changsha, Hunan 4 10081 ,P .R. China 


Abstract 

I have investigated the evolution of the field entropy in the two-photon JCM in the 
presence of the Stark shift and examined the effects of the dynamic Stark shift on the 
evolution of the field entropy and entanglement between the atom and field. My results have 
shown that the dynamic Stark shift plays an important role in the evolution of the field 
entropy in two-photon processes. 


1 Introduction 

The two-photon Jaynes-Cummings model [1] describing the interaction of a single-mode quantized 
field with a two-level atom through intermediate state involving the emission or absorption is one 
of the most intensively studied models in quantum optics . In this model, when the two atomic 
levels are couple with comparable strength to the intermediate relay level, the Stark shift becomes 
significant and cannot be ignored [2-5]. Puri and Bullough [3], A. Josh [4], Tahira Nasreen and 
Razmi [5] studied the influences of the Stark shift terms on the atomic inversion and dipole 
squeezing . Tahira Nasreen and Razmi [5,8] discussed the effect of the Stark shift on the Atomic 
emission and cavity field spectra in the two-photon JCM .These works have shown that the 
dynamic Stark shift plays an important role for the properties of dynamics in two-photon JCM . 
On the other hand, recently much attention has been focused on the properties of the entanglement 
between the field and atom in the Jaynes- Cummings model (JCM)[9-15].Phoenix and Knight [9] 
have shown that the partial entropy is a convenient and sensitive measure of entanglement between 
the atom and field . The time behavior of the field (atomic) entropy reflects time behavior 
of the degree of entanglement between the field and atom in JCM . The higher the entropy, 
the greater the entanglement, the information concerning the field is obtained by measurement 
performed on atoms .For the two-photon JCM, Phoenix and Knight [11], and Buzek [10] studied 
the evolution of the field entropy and the entanglement between the field and atom .The author 
[13] also examined the influence of atomic coherence on the evolution of field entropy in two- 
photon processes .However, these results are obtained in the case the Stark shift is ignored .In this 
paper, to make the two-photon JCM closer to the experimental realization,! include the effect of 
the dynamic Stark shift in studying the evolution of field entropy and entanglement . The results 
for the entropy evolution and entanglement incorporating the Stark shift are radically different 
from the results obtained in the absence of the Stark shift. 
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2 The reduced density operator and the field entropy 
calculation formalism for two-photon JCM in the pres- 
ence of the Stark shift 


In this paper, the model considered consists of a single-mode cavity field of frequency u> with 
an effective two-level atom of transition frequency u 0 through two-photon transitions in a lossless 
cavity. The excited and ground states of the atom will be designated by | +) and | — ), respectively 
. I assume these states to have identical parity, whereas the intermediate states, labeled | j ) 
0—3,4,...), are coupled to | +) and | — ) by a direct dipole transition and so located as to give rise 
to a significant Stark shift. The effective Hamiltonian describing such a model has form [3] 

He/f = u>a+a + u 0 S z + a + a{0 2 | +>(+ | +/? x | -)(- |) + g(a +2 S _ + a 2 S+), (1) 


where I have chosen units such that h — 1. a + and a are the creation and annihilation operator 
of the cavity field; S z =\ +){+ | - | -)(- |, £ + =| +)<- |, and =| -)(+ | are the atomic 

flopping operators. /5i and 0 2 are the parameters describing the dynamic Stark shift of the two 
levels due to the virtual transitions to the intermediate relay level, and g is the atom-field coupling 
constant . For simplicity ,1 consider on-resonance interaction ,so that ujq = 2u . By diagonalizing 
He// in the manifold of states | + , n) and | — , n + 2), the time-evolution operator in the interaction 
picture can be obtained [5] 


Ui(t) = 


U n in), U l2 (n) \ 
U 2 i(n), U 22 (n) J ’ 


( 2 ) 


where I have written 


Un(n) = sin 2 (0 n ) exp(— iA+i) + cos 2 (0„) exp(— iX^t) 

Ui 2 (n) = - sin(20 n )[exp(-zA+t) - exp(-iA^t)] = U 2 i(n) 

U 22 {n) = sin 2 (0 n )exp(-iA~t) + cos 2 (0 n )exp(-zA^t), (3) 

With 


sm (o n ) = i=(i + |r ) 1/2 


A± 


a 


9 


V2 

n(l + r 2 ) + 2 r 2 


2r 


±a 


= [5 2 (n+ l)(n + 2) + Q 1/2 

9 n_2l 


6* = ^[n(l — r 2 ) 


2r J 


r = (A/A) l/8 , (4) 

I consider the at time t=0 the atom is in a coherent superposition of the excited and ground states 


0 . V) = c °s(t|) I +) + exp(-zV) sin(^) | -) , 


( 5 ) 
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and the field is an arbitrary superposition of Fock states, so that at time t=0 the density operator 
for system 


OO OO Q Q 

p{ o) = EE F n ,m{ cos(-) I +,n) + exp(-z^) sin(-) | -,n)} 

n=0m=0 1 L 

9 6 

x (cos(-) (m, + I + exp(iyj) sin(-) (m, - |}. (6) 

where 9 is the degree of excitation, ip is the relative phase of the two atomic levels, F n<nx — F n F^ 
and F n are coefficients in the Fock- state. At any time t > 0 the reduced field density operator 
for the system is given by 

P/(t) = Tr a tom{0 J (t)p(0)Uj'(t)} 

1 OO OO 

= cos 2 (-60 £ E F n , m {Un(n)U* n (m) \ n)(m \ +U 2 i(n)U^{m) \ n + 2)(m + 2 |} 

^ n= 0 m= 0 


+ -sin(0)ex P M £ E F ntm {U u {n)U* n {m - 2) | n)(m - 2 


+U 2 i(n)U 22 (m — 2) | n + 2)(m |} 


+ -sin(0)exp(-i</j) E E F n ,Tn{U 12 (n - 2)f/* 1 (m) |n-2)(m 

^ n=0 m= 0 

+t/ 22 (n - 2)U 2x (m) | n)(m + 2 |} 


+ sin 2 (-0) E E - 2)t/* 2 (m - 2) | n)(m 

n=0 m = 0 

+[/ 21 (n-2)[/* 2 (m-2) \ n-2){m-2\}. 


The reduced density matrix Equation (7) has included the influences of the Stark shift and atomic 
coherence .Following the work of Phoenix and Knight [9], the eigenvalues and eigenstates of the 
reduced field density operator Equation (7) may be obtained 

*?(() = (tqt'i,) ± exp(=F«) I I 

= <f/ 2 f^ 1 >±exp(±6)|(t';^ 1 >l, (8) 

1^(0) = -i;"—— {expP(»g±g)] | C/i,)±exp[-f(t{±g)] | U^)} . (9) 

yJ2irf{t)cosh(S) 2 2 

where 

1 OO 1 OO 

| U' n ) = cos(-i 9) E F n U\i{n) \ n) + sin(-O) exp(-*V) E F n U l2 (n-2) | n - 2), 

^ n=0 ^ n=0 
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I ^ 21 ) = cos(;U) jr F n U 2X (n) I n + 2) + sin(^0) exp(-iv?) £ F n U 22 (n - 2) | rc)'. 

Z n=0 * n=0 

* = (10) 

6 = sinh~ 1 (x)- (11) 

I can obtained calculation formalism of the field entropy S/(t) in terms of the eigenvalue nf of 
the reduced field density operator 

S/(t) = -[nj(t)lnrf(t) + nj (t)lnnj (t)]. (12) 

The field entropy S/(t) given by Equation (12) shows that the field entropy depends not only on 
the field statistics parameters F n m ,but also on the Stark shift parameter r = (fi x / 0 2 ) x ^ 2 and the 
initial state of atom. Especially when 0 X = 0 2 = 0 ,Eq(12) give the results of Refs[l 1-13] for K=2. 

3 Numerical results 

I now discuss the numerical results for the field entropy S/(t) is given by Equation (17) when the 
initial field state is in the coherent state | a) 

F n , m = exp(— | a | 2 )tf*a"7(n!ro!) 1 ' 2 . (13) 

with a =| a | exp (ifi) and | a | 2 = n, where (3 is the initial phase of the field and n is the 
average number of photons in the coherent state. Here I hope to learn about the roles played by 
the Stark shift . The numerical results of equation (12) are show in Figs. 1-3 for different values 
of the Stark shift parameter r and different the initial states of the atom with n = 20 . 

; Figurel Figs.l Effects of the Stark shift on the evolution of the field entropy. 9 — 0 ,atom 
initially in excite state, field in the coherent state with mean photon numbers n = 20 ,(a),no Stark 
shift (0 X = 02 = 0 );(b),r = 1(0 X = 0 2 ) ;(c),r = 0.5 ; (d),r = 0.3 

; Figure2 Figs.2 Effects of the Stark shift on the 
evolution of the field entropy. 9 — 7t/2 , <£> — 2/3 = 0 ,atom initially in trapping state, the rest 
parameters as Figs.l. (a), no Stark shift (0 X = 0 2 = 0 ); (b),r = \(0 X = 0 2 ) ;(c),r = 0.5 ;(d),r = 
0.3 ; Figure3 Figs. 3 Effects of the Stark shift on the evolution 

of the field entropy. The same as Figs.2 but </? - 20 = 7t/ 2. (a), no Stark shift (0 X — 0 2 = 0 ); 
(b),r = l(/?i = p 2 ) ;(c),r = 0.5 ;(d),r = 0.3 3.1 Atom Initially in the excited state In 

Figs.l(a),I have the case 0 = 0 (i.e. the atom is in the excited state) and 0 X = 0 2 = 0(i.e. in 
the absence of the Stark shift), corresponding to the evolution of the field entropy in the standard 
two-photon JCM obtained by Ref[l 1-1 3] .1 note that the field entropy evolves at periods 7 r/< 7 ,when 
t = 7i7 r/< 7 , (n = 0, 1, 2, 3...),S/(t ) evolves to the zero values and the field is completely disentangled 
with the atom, while when t = (n + l/2)ir/g. Sf(t ) evolves to the maximum value, and the field is 
strongly entangled with the atom . The results for the evolution of the field entropy S/(t) in the 
presence of the Stark shift are plotted in Figure l(b)-(d). In Figs. 1(b), the Stark shift parameter r 
is given as 1 (namely 0 X = 0 2 ), this correspond to the case the two levels of the atom are equally 
strongly couple with the intermediate relay level . By making a comparison between Figs. 1(a) 
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and Figs. 1(b), I find that the evolution of the entropy is almost similar for both cases. This 
result corresponds with the fact that in two-photon processes, the Stark shift creates an effective 
intensity dependent detuning A N = p 2 - Pi [16]. When r = 1.0 , (namely Pi = /? 2 )thus An = 0.0 , 
the Stark shift does not affect the time evolution of the field entropy . In Figs. 1(c), I show the 
case r=0.5,in which the two levels have the unequal Stark shifts (Pi < P 2 ). I note that the Stark 
shift leads to decreasing the values of maximum field entropy and increasing the values of field 
minimum entropy .It also results in increasing frequency of the field entropy vibration . As the 
parameter r further decreased (e.g.,r=0.3,see Figs. 1(d)), the values of the maximum field entropy 
and the degree of entanglement of the field-atom further reduced . 3.2 Atom initially in the 
superposition states When 6 = 7r/2 and <p — 2/? = 0 , the atom is initially in trapping state, 
the evolution of the field entropy are plotted in Figs. 2 . Figs.2(a) show the evolution of the field 
entropy in the absence of the Stark shift while Figs. 2(b) show the case that r=1.0 .1 can see that 
under the condition the atom is initially in trapping state, the field entropy obviously reduced with 
comparison Figs.l (a),(b) in where the atom is initially in excited state, and the evolution of field 
entropy in the case r=1.0 is almost same as that in the absence of the Stark shift .Figs.2(c) and 
Figs.2(d) show the evolution of Sf(t) in two cases that r=0.5 and r=0. 3, respectively. As is visible 
from the figures, the effects of the dynamic Stark shift are more pronounced when r deviates 
from unity .On the other hand .when atom in initially trapping state, as r decreased, the values of 
maximum field entropy increased, indicating that the Stark shift leads to increasing the degree 
of entanglement between the field and atom, which is contrary to the case the atom is initially 
excited state . The results for the evolution of the field entropy S/(t) as 6 = 7t/2 , ip — 2/3 = 7r/2 
and various values of Stark shift parameter r are presented in Figure 3 .In these cases, I note 
that the maximum field entropy always remains at its maximum values, regardless of the chosen 
value of r . A possible explanation for above behavior of field entropy evolution can be performed 
in terms of the Bloch vector in semi classical theory [17] in the next section . 

4 Semiclassical interpretation of the evolution behavior 
of the field entropy 

I can find that the larger the extent of the Bloch vector’s motion,the greater the values of the 
maximum field entropy by examining the evolution of field entropy in semiclassical version . By 
replacing the field annihilation operator a by the c number v = \/Eexp{i(ut— p)\ I get semiclassical 
version of eq.(l) 


H e ff — nu + UoS z + n(P 2 | +)(+ I +pl | —){— |) 

+ ng(S + exp[— i2(ut — /?)] + SL exp[i2(ut — /?)]) . 


The motion equations of operators = +, — , z)can be written 


dS+ 

dt 


= iu>oS + — i2gn exp [i2(ut — P)]S Z + inANS + 


> 


dS_ 

dt 


— iuJoS _ + i2gnexp[—i2(ut — P)\S Z — inA^5_ , 


(14) 


(15) 

(16) 


411 


( 17 ) 


dS , 
d£ 


= ign(exp[i2(ut — 0)}S- — exp[— i2{ut — /?)]S + ) , 


where 




02 ~ 01 = 


5(1 - r 2 ) 


(18) 


is an effective detuning created by the Stark shift . Define the two slowly varying operators which 
involve the coherence between the two atomic states 


S x = ^{exp[-z(a;o£ - $)]£+ + exp[i(wo£ - $)]£_} , (19) 

S y = ^{exp[-z(o;oi - ^ > )]5'+ - exp(i(wbf - . (20) 

2x 

where <F is a phase angle that may be chosen at will . The atom interacting with the field obeys 
the optical Bloch equation 

^ = n„(() x s(t). ( 21 ) 

Where S(t) is the Bloch vector for the atom 

S(t) = {S,(t),S y (t),S,m ■ (22) 


and Cl s (t) is the driving field vector, which can be written by using eqs.(14)-(22) 

fi a (t) = {2gncos[(2u; — uj 0 )t — ( 20 — 4>)], 2gnsin[(2u — Uo)t — ( 20 — 4>)],nA^} . (23) 

Generally S(£) precesses in a cone about fi s (t). The extent of Bloch vector’s motion is largest when 
S(t) and S"2 s (£) are orthogonal and minimum when S(£) and f2 s (£) axe parallel or antiparallel. Thus 
the time evolution of the Bloch vector is quite different for different initial preparations. If the 
atom is initially coherent superposition state given by eq(5), the Bloch vector at time t=0 can be 
expressed [17] 

S(0) = sin(0) cos(4> - ip), ^ sin(0) sin($ - ip), ^ cos(0 )} , (24) 

Cl s (0) = {2f/ncos(4> — 20),2gfism($ — 2/?),nA^} . (25) 

For simplicity, I let 4> = <p,the initial Bloch vector in this case is on the x-z plane. When the 
atom is initially in excite state (5 = 0 ),the initial Bloch vector S(0) = {0,0, and the vector 
n s (0) = {2gn, 0, nA;v} . Under this condition, the extent of Bloach vector’s motion is dependent 
of the effective detuning A/v created by the Stark shift . When 0 { = 0 2 = 0 (in the absence 
of the Stark shift) or r = 1.0(/?i = 0 2 ) , A n = 0 , the vector S(0) and fl s (0) are orthogonal 
and the extent of Bloach vector is the maximum for these two cases so that the maximum field 
entropy always remains at its maximum value and the degree of the entanglement of atom-field is 
the largest (see Figs. 1(a) and (b)). With r reduced, Aw increased, S(0) and fi s (0) are no longer 
orthogonal and the extent of Bloach vector’s motion is reduced .This leads to the values of the 
maximum field entropy and the degree of entanglement of atom-field decreased as the Stark shift 
parameter r reduced under the condition the atom is initially in excite state . When the atom 
is initially trapping state (namely 0 = — 20 = 0), Bloach vector S(0) = {i, 0, 0}, f2 s (0) = 

{2gn, 0,hAw}.The motion extent of S (t) is also dependent of the effective detuning A^r . When 
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fa = fa — Q (in the absence of the Stark shift) or r = 1.0(/?i = fa) » An = 0-0 an d the vector S(0) 
is parallel to f2 s (0),the motion extent of the S (t) tends to zero, and the values of the maximum field 
entropy obviously reduced (see Figs. 2(a), (b). With parameter r decreased, A N increased and the 
vector S(0) and fi s (0) are no longer parallel . The greater the effective detuning An (the smaller 
the parameter r), the larger the extent of the S(f) ,this leads to the values of maximum field 
entropy and the degree of entanglement of the atom-field increased as parameter r decreased (see 
Figs. 2(c), (d)) under the atom is initially in trapping state . Furthermore, when 9 = n/2 Relative 
phase </? — 2/? = 7 t /2 is chosen, S(0) = {£, 0, 0}, O s (0) = {0,2gn,nA N } are completely orthogonal 
irrespective of values of An . Therefore, at this case, the extent of motion and the values of 
maximum field entropy are independent of the influences of the Stark shift, and are maximum at 
all(see Figs.3(a)-(d)). 
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The figure captions 

Figs.l Effects of the Stark shift on the evolution of the field entropy. 0 = 0 ,atom initially 
in excite state, field in the coherent state with mean photon numbers n = 20 ,(a),no Stark shift 
(ft = ft = 0 );(b),r = l(ft = ft) ;(c),r = 0.5 ;(d),r = 0.3 

Figs.2 Effects of the Stark shift on the evolution of the field entropy. $ = ir/2 , <p — = 

0 ,atom initially in trapping state, the rest parameters as Figs.i. (a), no Stark shift (ft = ft = 0 ); 
(b),r = l(ft = ft) ;(c),r = 0.5 ;(d),r = 0.3 

Figs.3 Effects of the Stark shift on the evolution of the field entropy. The same as Figs.2 
but — 20 = 7r/2. (a), no Stark shift (ft = ft = 0 ); (b),r = l(ft = ft) ;(c),r = 0.5 ;(d),r = 0.3 
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